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1. Introduction 



The basic definition of a fusion coefficient is that it gives the number of independent 
couplings between three different fields in conformal field theory (cf. also the introduction 
of [1]; for a review of conformal field theory, and in particular fusion rules, see [2]). Even in 
theories with a Lie group symmetry, the so-called Wess-Zumino-Witten (WZW) models, 
an intrinsic conformal-field theoretical characterisation is unavoidable. This is manifest in 
formulae for the fusion coefficients: the most fundamental one is the Verlinde formula [3] , 
that expresses a fusion coefficient in terms of modular S matrix elements: 

c c c* 

Here we use notation appropriate to a WZW model in which primary fields are in one-to- 
one correspondence with the integrable representations of the spectrum-generating affine 
algebra at a fixed level k (this set is denoted by P£) and stands for the basic represen- 
tation, whose finite projection is the scalar representation. Fields are not distinguished 
from their representation labels. The matrix S specifies the linear modular transformation 
properties of the characters of the primary fields among themselves. Up to a constant fixed 
by unitarity, it takes the form 

- ] -—(w(X + p)^ + p)) (1.2) 

where g stands for the dual Coxeter number of the algebra under consideration, p is the 
Weyl vector, A is the finite projection of the affine weight A and W is the finite Weyl group. 

The remarkable fact that the ratio of two S matrix elements is a finite character 
evaluated at a special point yields a close relation between fusion and tensor-product 
coefficients. Indeed, since the finite character and its evaluation read 

XA= E„ w «»)<»' and xx{i)= E^W^ < L3 > 

we observe that 

Xx(Z) = ^ L with Z = --^( a + p) (1.4) 
£>o,a k + g 

This leads to the Kac- Walton formula which relates the fusion and the tensor-product 

coefficients. 
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The Verlinde formula does not make manifest the basic integrality property of the 
fusion coefficients. The S matrix elements being in general complex numbers, it is not 
even clear at first sight that the fusion coefficients are real (this follows from the unitarity 
property of S) . The integrality is ensured by the Kac- Walton formula, but in this case the 
positivity is not manifest. 

It is mainly with the aim of displaying manifestly non-negative formulae for fusion 
rules that we have looked for fusion generating functions [4]. Although the construction 
of explicit generating functions has an intrinsic interest, we regard the unravelling of the 
concept of threshold level - reviewed below - as being the most important outcome of 
this analysis. It leads to a complete characterisation of fusion coefficients in terms of the 
corresponding tensor-product coefficients and a set of threshold levels. 

As a result, the interest has shifted from the construction of fusion generating functions 
to the search for threshold-level computing techniques. For su(N), N = 2,3,4, it has 
been found that the threshold level is coded in a simple way in the Berenstein-Zelevinsky 
triangles [5] (cf. also section 7.1 of [1]) describing the various distinct couplings of a tensor 
product [6, 7]. However, these formulae are difficult to generalise to larger values of N. 
Moreover, this approach, based on a diagrammatic description of the tensor product, is 
limited to the su(N) algebras. 

The aim of the present paper is to apply the machinery developed in [1] to these 
problems. We find new generating functions for sp(4) and sw(4), together with a closed 
form expression for their threshold levels. More importantly, we introduce the concept 
of fusion basis, that is, the set of linear and homogeneous Diophantine inequalities that 
describes completely the fusion rules. 

The article is organised as follows. In section 2, after introducing some notation, we 
present a brief review of fusion rules and show, with the example of su(2), how tensor- 
product generating functions and the Kac- Walton algorithm can be used to construct 
fusion-rule generating functions. A more powerful approach to the problem is then elabo- 
rated in section 3. It relies on the conjectural existence of a linear and homogeneous set 
of inequalities that provides a complete description of fusion rules. Given a set of fusion 
elementary couplings, Farkas' lemma is then used as a technique to extract the underlying 
inequalities. This is what we call a fusion basis, i.e., the basis in terms of which these fusion 
elementary couplings are the elementary solutions. A complete analysis of the sw(3), sp(4) 
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and su(4) cases is presented in section 4, 5 and 6 respectively. In all three cases, the 
general expression for the threshold levels is obtained explicitly. Various arguments (based 
on Giambelli-type formula and level-rank duality) supporting our results are presented in 
Appendix A. In Appendix B, we recall previous conjectures and clarify their relation to 
those formulated here. 

2. Fusion rules 

Let g be the affine Lie algebra corresponding to the finite Lie algebra g. Quantities with 
hats generally refers to g. The fundamental weights of g are denoted by u)i, % = 0, 1, r, 
where r is the rank of g. An affine weight may be written as 

r 

A = J]A i a; i = [Ao,Ai,...,A r ] (2.1) 

i=0 

If the Dynkin labels Aj are nonnegative, the weight A is the highest weight of an integrable 
representation of g at level k, with k defined by 

r 

k = Y. X ^ ( 2 - 2 ) 

The a( are the co-marks: o,q = 1, and the remaining a( are the coefficients of expansion 
of the longest root of g in terms of the simple coroots. The set of such weights is denoted 

To the affine weight A, we associate a weight A of the finite algebra g 

r 

A = y^^jUj = (Ai,...,A r ) (2.3) 

where u>i for (i = 1, r) are the fundamental weights of g. A is thus uniquely fixed from 
A and k. The set of integrable finite weights is written P + . 

In the conformal field-theory context, fusion rules yield the number of independent 
couplings between three given primary fields. Here we are interested in fusion rules in 
WZW models [8,9], whose generating spectrum algebra is an affine Lie algebra at integer 
level. 
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Denote the multiplicity of the representation v in the fusion rule A x ft by 



A x £ = £ A/$ * * (2.4) 



and denote by Nx^ the multiplicity of the representation v in the tensor product A ® fx: 



A <g> n = ^ AV 1/ (2.5) 



where by abuse of notation, we use the same symbol for the highest weight and the highest- 
weight representation. The precise relation between tensor-product and fusion-rule coeffi- 
cients is given by the Kac-Walton formula [10,11,12]: 

^ x f= E A/yu«o (2.6) 

w is an element of the affine Weyl group W, of sign e(w), and the dot indicates the shifted 
action, 

r 

w ■ A = w(X + p) — p p = ^^Qi (2-7) 

The Kac-Walton formula can be transformed into a simple algorithm: one first calcu- 
lates the tensor product of the corresponding finite weights and then extends every weight 
to its affine version at the appropriate value of k and shift-reflects back to the integrable 
affine sector those weights which have negative zeroth Dynkin label. Weights that cannot 
be shift-reflected in the integrable sector are ignored (for example this is the case for those 
which have zeroth Dynkin label equal to —1). 

The affine extension of the weights that occur in the tensor product may not be inte- 
grable at level k but are integrable at level 2k. If we divide the weight space into domains 
that are mapped into each other by the application of the affine Weyl reflections, then the 
affine reflections which contribute to the Kac-Walton algorithm, apart from the identity, 
are those corresponding to the domains next to the fundamental alcove and which lies in 
the P + cone. This is a crucial property of the Kac-Walton algorithm for its application 
to the construction of fusion-rule generating functions. Let us denote by Wf this finite 



subset of the affine Weyl group that need to be considered . For instance, the elements of 
Wf for the lowest rank algebras are: 



su(2) : 


W f 


= {id, s } 


su(3) : 


W f 


= {id, s , s\Sq, s 2 sq} 


Ht(4) : 


w f 


= {id, s , stso, s 3 s , S2S1S0, S2S3S0, S1S3S0, s sis 3 s } 


sp(4) : 


Wf 


= {id, s , sis , s sis } 


G 2 : 


Wf 


= {id, So, SxSq, S 2 S , SoS2S!S } 



where Sj denotes the reflection with respect to the root ctj. This set of elements w can be 
characterised as follows: these are the elements w of the affine Weyl group that satisfy the 
requirement: 

w{2a% + a\ + ■ ■ ■ + 0%, a\, ■ ■ ■ , a y r } E (2.9) 

where stands for the set of positive real coroots of the affine algebra under consideration 
and r stands for its rank. This condition is adapted from [13] as further analysed in [14]. 

Note also that (2.6) may be rewritten as: 

v G : J^°= E W<*>) (2-10) 

where it is understood that w ■ v stands for its finite part since it is an index of the tensor- 
product coefficient. This allows us to study in isolation the contribution of a single weight 
in the fusion. For instance, for su(2) that reads 

^ *=V- S °' P = A/"a iMi Ul ~ -A/a iMi 2k+2 ~ Ul (2.11) 

Here is an illustrative example of the Kac- Walton algorithm that will also serve to 
introduce the key notion of threshold level. Take the following sp(4) tensor product: 
(1,1)®(1,1). Its decomposition reads 

(1, 1) ® (1, 1) = (0, 0) © (0, 1) © 2 (2, 0) © (0, 2) © (0, 3) © 2(2, 1) © (2, 2) © (4, 0) (2.12) 

The sp{A) comarks are all equal to one so that the affine extension of a weight (m, n) at level 
k is [k—m—n, m, n]. At level 2, the weights (0, 3) and (2, 1) are ignored (they have v Q = — 1) 
and the remaining non-integrable weights are [—2,2,2] and [—2,4,0]. Since the zeroth 
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simple root is a = [2, -2, 0], we have s ■ [-2, 2, 2] = [0, 0, 2] and s ■ [-2, 4, 0] = [0, 2, 0], 
so that the resulting fusion is 



[0, 1, 1] x [0, 1, 1] = [2, 0, 0] © [1, 0, 1] © [0, 2, 0] (2.13) 

In the above example, we see that the weights (0, 0), (0, 1), (2, 0) appear first at level 
2. It is easily checked that they reappear at every level k > 2. We then say that their 
threshold level, usually denoted by ko, is 2. The threshold level is thus the smallest value 
of k such that the fusion coefficient J\f~ k } v is non-zero. If we indicate the threshold level 
by a subindex, by considering the extension of the above tensor product at different levels, 
we find 

(1, 1) © (1, 1) = (0, 0) 2 © (0, 1) 2 © (2, 0)2 © (2, 0)3 © (0, 2)3 

(2.14) 

© (0,3)3 © 2(2,1)3 © (2,2)4 © (4,0)4 

To read off a fusion at fixed level k, we only keep terms with index not greater than k. The 
concept of threshold level was first introduced in [4]. It can be shown (cf. ref. [6]) that 
the existence of a threshold level is a consequence the depth rule of Gepner and Witten 
[9] . The notion of threshold level implies directly that 

M k) ° < M k+1) ° and lim J&) p = JV A /. (2.15) 

To the triplet (A, /i, v) there corresponds distinct couplings, hence values of 
ko, one for each distinct coupling. Let us denote these by k$\i = 1, Nx^ , implementing 
in this notation the natural ordering k^ < k^ +1 \ Then 



. max(z) if k > k {l) and A/" A / ^ 
M } " = \ " m (2-16) 

1 if £;<£;o (1) or = 0. 



Further variations on the idea of threshold level are presented in [15]. 

Let us finally note that the fusion coefficients are invariant under the following action 
of the outer-automorphism group [10] 



<V AA ' P =^° (2-17) 

AX,A'n Xfi v ' 
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For example, for sp(4), the non-trivial outer automorphism a exchanges the zeroth 
and second root, or equivalently, it acts on weights as a[Ao,Ai,A2] = [A2,Ai,Ao]. Acting 
on the fusion (2.13) as 

a[0, 1, 1] x a[0, 1, 1] = [1, 1, 0] x [1, 1, 0] = [2, 0, 0] © [1, 0, 1] © [0, 2, 0] (2.18) 

which is easily checked from the tensor product 

(1, 0) © (1, 0) = (0, 0) © (0, 1) © (2, 0) (2.19) 

which is non-truncated at level 2. Other fusions at level 2 can be obtained from (2.13) by 
acting on the weights as follows 

a[0, 1, 1] x [0, 1, 1] = [1, 1, 0] x [0, 1, 1] = a[2, 0, 0] © a[l, 0, 1] © a[0, 2, 0] 

(2.20) 

= [0,0,2] © [1,0,1] © [0,2,0] 

The algorithm underlying the Kac- Walton formula suggests a simple road to the con- 
struction of fusion-rule generating functions, that is by starting from the tensor-product 
calculation, but keeping track of the level and taking into account the action of the affine 
Weyl group. We illustrate the method for the simple su(2) case. Recall that the su(2) 
tensor-product generating function reads 



QSU{2) M ^ = 1 

V ' ' ; (1- LM)(1- LN)(1- MN) y 1 

We start with the generating function 

F (ti ,L.M.iV) = (i _ tQ(i _ LM)( ;_ £Ar)(i _ MAr) . (2.22) 

This is just the generating function for su(2) tensor products divided by (1 — d). The 
exponent of d will be identified with the level. We will proceed to the generating function 
for su{2) fusion rules by modifying (2.22). First note that at level k we need only consider 
the products of su(2) representations (a) with a < k. The generating function (2.22) 
includes products of representations which violate this condition. To keep terms of the 
form d k L a with a < k introduce a dummy variable x (using MacMahon's notation - cf. 

Q - —-F(dx 7 Lx' 1 , M, N) (2.23) 

= (1-x L ) 



This first converts d k L a to x~ m+k ~ a d k L a , with m > and then keeps terms of degree 
zero in x which corresponds to keeping the terms of F(d, L, M, N) with a < k as required. 
This yields: 

1 (2.24) 



(1 - d)(l - dLM)(l - dLN)(l - MN) ' 
Repeating this procedure with L replaced by M yields: 

1 — d 2 LMN 2 

G(d, L, M, N) = (1 _ _ d LM)(l - dLN)(l - dMN)(l - dLMN 2 ) ' (2 ' 25) 

This is still a generating function for tensor products, but with the size of the representation 
Dynkin labels restricted to be less than or equal to the level. 

To take into account the affine Weyl group, consider a term in the expansion of the 
generating function which contains d k N c . If c > k + 1 then this representation is reflected 
back into the fundamental region of the affine Weyl group: c i— > c— 2(c — k — 1) = — c+2k+2 
or d k N c i— > d k N 2k ~ c+2 . Since this is a reflection, the corresponding character must be 
subtracted. In principle other affine Weyl transformations might be necessary to obtain 
a weight in the fundamental domain, but, as discussed earlier, for su(2) one reflection 
suffices. At the level of generating functions the effect is to replace 

G(d, L, M, N) ^ G(d, L, M, N) — N 2 G(dN 2 , L, M, iV -1 ) (2.26) 

Note that the new generating function contains terms with negative powers of iV and also 
terms with c > k. To obtain the final function we projected out the required terms as 
above. Although this calculation is somewhat long (the verification here was done on a 
computer), the final result is very simple: 

G ?u{2) = 1 , 22 7) 

(l-d)(l-dLM)(l-dLN)(l-dMN) y 1 

This has first been written down in [4]. There are thus four elementary couplings: 
E :d : (0) ® (0) D (0)i E 2 : dLN : (1) ® (0) D 

(z.zo J 

Ei : dLM : (1) <g> (1) D (0)i, E 3 : dMN : (0) <g> (1) D (l)i. 

As explained above, subscripts indicate the threshold level. 

The notion of a model was discussed in [1]. A model for this generating function is 
Q[E , Ei, E 2 , E 3 ) with the gradings of E , • ■ • , E 3 respectively given by (1, 0, 0, 0), (1, 1, 1, 0), (1, 1, 0, 1) 
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and (1, 0, 1, 1) for the ordering X = d,X 1 = L, X 2 = M, X 3 = N. As for the finite su(2) 
case, there are no relations between the elementary couplings. 

The generalisation of the above calculation to other affine Lie algebras is straight- 
forward. Starting from the tensor-product generating function augmented by the factor 
1/(1 — d), where d keeps track of the level, one first enforces the integrability requirement 
of the first two weights (those that are fused together) ; one then implements all the affine 
reflections of the set Wf on the third weight and projects the alternating sum onto the 
integrable sector. However, even though the strategy is clear, the computations become 
rapidly very complicated. 

To bypass this difficulty, we have argued (cf. [1] section 3) that the use of a direct 
description of tensor products in terms of a system of inequalities (e.g., the Littlewood- 
Richardson (LR) inequalities underlying their combinatorial description for calculating 
tensor products - cf. [1] section 4) simplifies the general procedure to a very large extent 
in addition to allowing us to use powerful algebraic results. We now look for a similar 
procedure here. However, this program faces an immediate difficulty since even for su(N), 
a combinatorial description of fusion rules is not known. Our method is instead to find an 
independent route leading to the elementary couplings. Indeed, the elementary couplings 
are really what we need in order to apply our Grobner basis machinery. Quite remarkably, 
it turns out that once elementary couplings are found, there is a method that allows us to 
reconstruct the underlying system of Diophantine inequalities. 

3. Fusion-rule elementary couplings 

The construction of this section depends upon the following: 

Fundamental conjecture: There exists a fusion basis, that is, a set linear and homogeneous 
inequalities involving k and containing as a subset, a tensor-product basis. 

For instance, the LR basis is a set linear and homogeneous inequalities. Every solution 
can be expanded in terms of the elementary solutions of these inequalities. For su(N), 
the conjecture amounts to the existence of a set of additional inequalities involving the 
level k that provide the proper truncation describing the fusion rules. The relation of this 
conjecture to the conjectures presented in [4] is discussed in the Appendix B. 
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Note that homogeneity is the key property which allows us to reconstruct the fusion 
basis from a set of fusion elementary couplings using Farkas' Lemma. This condition does 
not necessarily hold, for example we have found that the Lie superalgebra osp(l,2) does 
not have a homogeneous basis. 

Given homogeneity and Farkas' lemma, the problem is reduced to finding a set of 
fusion elementary couplings. The Kac- Walton algorithm is one possible approach, but a 
rather difficult one. Instead, we will introduce a simpler approach based on the outer- 
automorphism group. Unfortunately, it relies on another conjecture. 

Let us start from the set of tensor-product elementary couplings {Ei, i £ 1} for some 
set / fixed by the algebra under study. For each Ei, we calculate the threshold level 
ko(Ei). This information specifies the affine extension of Ei. The affine extension of a 
tensor-product elementary coupling is necessarily a fusion-rule elementary coupling given 
our hypothesis that the fusion basis contains, as a subsystem, the set of inequalities that 
describe tensor products. Denoting by a hat the affine extension of a tensor-product 
elementary coupling 

Ei = d ko ^Ei (3.1) 

we have then a partial set of fusion elementary couplings with the set {Ei,i e /}. Our 
conjecture is that the missing fusion elementary couplings can all be generated by the 
action of the outer-automorphism group whenever this group is nontrivial: 

The outer- automorphism completeness conjecture: The complete set of elementary cou- 
plings {Ei, i G J} for a set J D I can be generated by the action of the outer-automorphism 
group on the set {Ei, i e /}, i.e., the full set is contained in {AEi}: 

{E t ,ieJ}c{AE t ,iel} (3.2) 

The action of the outer-automorphism group on a coupling is defined as follows. Let 
the three weights in the coupling be {A, fi; u} where z> C A x fi, then 

A{\, fi; u} = {AX, A' ft; AA'u} (3.3) 

where A, A' are arbitrary elements of the outer-automorphism group; the conjectured 
completeness requires the consideration of all possible pairs (A, A') . 
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It should be stressed that we do not suppose that the action of A on an elementary 
coupling will necessarily produce another elementary coupling. Indeed, the resulting cou- 
pling could be a product of elementary couplings. What is conjectured here is that all 
fusion elementary couplings can be generated in this way. 

If the outer-automorphism group is trivial, we expect that there will a single extra 
elementary coupling, the one associated to the scalar coupling: E . 

As a simple example consider su(2). Start with the elementary coupling E\ : (1) ® 
(1) D (0). It is easy to show that this coupling arises at level 1. This is thus the value 
of its threshold level. The corresponding fusion is [0, 1] x [0, 1] D [1, 0]. We now consider 
all possible actions of the outer-automorphims group on it. Since this group is of order 2, 
there are 4 possible choices for the pair 



with a[Ao, Ai] = [Ai,Ao]. This generates the following set of four elementary couplings 
found previously (cf. eq (2.28)): 



Let us then suppose that we have a complete set of fusion elementary couplings which 
are the elementary solutions of set of linear and homogeneous inequalities that we are 
looking for. A standard theorem in the theory of linear Diophantine equations (cf. [16]) 
states that every non-negative integer solution of a given set of homogeneous Diophantine 
inequalities for the variables Xi (e.g., for su(N), these are the {Xi,riij}) can be generated 
from a non-negative combination of the fundamental solutions. Hence, given the set of 
elementary couplings any coupling can be decomposed (maybe not uniquely) in the 

form Yli^P- Let the grading variables representing the Xi be denoted by X^. To the 
expression of g{Ei) corresponds a vector of components e^ , which is the vector form of 
the elementary solutions of the Diophantine equations. In other words, 



(A,A')e{(a,a), (1,1), (l,a),(a,l)} 



(3.4) 



E :d : [1,0] x [1,0] D [1,0] 
E x : dLM : [0,1] x [0,1] D [1,0] 



E 2 : dLN : [0,1] x [1,0] D [0,1] 
E 3 : dMN : [1,0] x [0,1] D [0,1]. 



(3.5) 



%: g(E i ) = Hx € j 



(3.6) 



3 
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Reading off a particular coupling means that we are interested in a specific set of non- 
negative integers {x{\ given by 



in terms of no n- negative integers c^. We are thus looking for the existence conditions for 
such a coupling. This is related to Farkas' lemma [17,18]. The standard, rational, form of 
the lemma is (cf. [17], corollary 7. Id): 

Farkas' lemma: Let V an m x n matrix with rational entries and let x G Q m . Then there 
exists a > 0, a G Q n such that Va = x if and only if for all u G Q m , u T V > implies 
u T x > 0. 

We can relate this to our problem in the following way. First note that the condition 
that for all u G Q m , u T V > implies u T x > is equivalent to the condition that for 
all u G Z m , u^V > implies u T x > 0. Necessity is clear and sufficiency follows since if 
u G Q m and u T V > then u = cu' with c G Q, c> and u' G Z m . Then u' T V > 0, so 
-u /T x > and multiplying by c gives the required inequality. 

Now consider the inequalities 



By writing u\ = Wi — Vi, Wi,Vi G N, i = 1 . . . m, we obtain a new system of linear Dio- 
phantine inequalities. It is not difficult to see that every solution to (3.8) can be obtained 
from a solution to this new system. Moreover, the new system of linear Diophantine in- 
equalities has a finite set of fundamental solutions. These give rise to a set of fundamental 
solutions to (3.8) such that every solution to (3.8) is a linear combination of these funda- 
mental solutions with non-negative integer coefficients. Call these fundamental solutions 
Si, i = l...k. Thus the condition that for all u G Z m , u T V > implies u T x > is 
equivalent to the condition sjx > 0, i = 1 . . . k. 

Putting all this together we obtain the following variation of Farkas' lemma: 

Lemma: Let V be an m x n matrix with rational entries and let x G Q m . Then there exists 
a > 0, a G Q n such that Va = x if and only if sjx > 0, i = 1 . . . k where s$, i = 1 . . . k are 
a fundamental set of solutions of the system u T V > 0, u G Z m . 

We can reformulate this Lemma over the integers in a form which is more convenient 
for our application: 




(3.7) 



u 1 V > 0, ueZ 



m 



(3.8) 
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Proposition A: Suppose V G M m;n (N) and let x G N m , a G N n . Then Va = x if and only 
if = ct^a, * = where itj, aij, z = 1 . . . are a fundamental set of solutions of the 

system u T V = a T , u G Z m , a G N n . 

To show this, suppose that Va = x and that u T V = a T with u G Z m , a G N n . Then 
u 1 x = u T Va = a T a. In particular this is true for the fundamental solutions. 

Conversely, suppose ujx = aja for every fundamental solution. Then u T x = a T a 
for every u G Z m , a G N n such that u T V = a T . Since V G M m>n (N), one set of solutions 
to u T V = a T , u G Z m , a G N n is given by taking u to be a suitable unit vector and a T to 
be a row of V which gives Va = x as required. 

To link the lemma to the situation presented above, we note that the entries Vij 
of the matrix V are given here by the numbers eji appearing in (3.6). Our analogue of 
the relation V a = x describes a generic coupling and our goal is to find the defining 
system of inequalities underlying the existence of this coupling. The equalities ujx = aj a 
i = 1 . . . k imply that x satisfies ujx > i = 1 . . . k since and a are non-negative. In 
general these inequalities have solutions which are not solutions of the former equalities 
for any a. For example if V = (2), then Va = x is 2a = x which is also the equality 
obtained from the second part of the Proposition A. Thus a; is a non-negative even integer. 
But the corresponding inequality is x > 0. However, we have found that for the particular 
systems we consider, this does not happen - as can be easily verified by computing the 
fundamental set of solutions to the inequalities ujx > i = 1 . . . k and verifying that they 
are the columns of V. 

As a simple illustration of this construction, let us work out the example of 
su(2). We use the LR variables {k, Ai, nn, n^} and the corresponding grading variables 
{d, Li, Nn, Ni 2 } in terms of which the elementary couplings and the corresponding vectors 
are 

E :d e = (1,0,0,0) 

E^.dL^ ei = (1,1,0,1) 

E 2 :dL x e 2 = (1,1,0,0) 

E 3 :dN 11 e 3 = (1,0,1,0) 

For future reference, we display the LR tableaux of the corresponding tensor-product 
elementary couplings 

Ei: T , E 2 : □ , E 3 : [T] (3.10) 
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To the fusion elementary couplings, we associate the vectors ej which form the matrix V 
with components Vij = e^: 

/I 1 1 IX 



V = 



110 
1 

Vo i o o/ 



(3.11) 



and so we have the matrix equation 

Va = x (3.12) 

This equation describes a general fusion coupling. We now want to unravel the underlying 
system of inequalities. For this, we use Proposition A, i.e., we find the fundamental 
solutions of u T V > 0. This is first transformed into a set of equalities u T V = a T by 
introducing new non-negative parameters af. 

Uq = «o Uq + U\ = OL 2 

(3.13) 

U + Ui+U 3 = «i Uq + U 2 = «3 

We next apply the vector-basis arguments (see Section 7 of [1]). Let us choose the as 
our independent variables. (This example is somewhat misleading due to its simplicity: in 
general not all the can be taken as the independent variables.) The dependent variables 
read then 

u = a u 2 = a 3 — a 

(3.14) 

U\ = OLl — OCQ U2, = CX\ — «2 

The 4 basis vectors are obtained by setting successively one ctj equal to 1 and all the others 
equal to 0. These vectors are written as and their entries are 

ei = (uo(ati = 1), ui(ati = 1), -u 2 («i = 1), u 3 (ai = 1); a , ati, a 2 , a 3 ) (3.15) 

With % = 0, 1,2,3, we find 

e = (1, -1, -1, 0; 1, 0, 0, 0) e 3 = (0, 0, 1, 0; 0, 0, 0, 1) 

d = (0, 0, 0, 1; 0, 1, 0, 0) e 2 = (0, 1, 0, -1; 0, 0, 1, 0) 



(3.16) 



These ei are manifestly linearly independent and they are non-negative expressions in the 
aj. In other words, their grading re-transcription of the above vectors (with Ui and A% 
denoting the grading variables of itj and ctj respectively) reads 

S = UoU^U^Ao S 2 = U 1 U 3 1 A 2 

£ 1 = U 3 A 1 S 3 = U 2 A 3 
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Here we see that all Si contain positive powers of the Ai (this is not generic and it reflects 
the simplicity of the su(2) case). Hence, all solutions are generated freely from the non- 
negative powers of the £j. 

The corresponding linear system of Proposition A is ei(x, — a) T = with x = 
(k, Ai, nn, "-12) and a = (a\, 0,2, 03, (14) non-negative integers: 

k — Ai - nn = 0-1 Ai - n\2 = «3 

(3.18) 

n i2 = a 2 nn = «4 

which are equivalent to the inequalities: 

k > Ai + nn Ai > n i2 

(3.19) 

ni2 > nn > 

The last three conditions define the LR basis. The first one is the additional fusion con- 
straint. 

In general, we will work the elementary solutions in their exponential version Si to 
keep the notation more compact and it should be clear that the (in) equalities can be read 
off as easily at this level. 

The construction of the su(2) generating function is now straightforward: since there 
are no relations between the elementary couplings, the generating function is simply (2.27), 
that is 

3 1 

G su{2) = "Q ^_ (320 ) 

i=o (! - E i) 

From the /c-inequality of the su(2) fusion basis, we read off the threshold level of a 
coupling as ko = Xi + nn, that is 

ko = (\ 1 + m+v 1 )/2 (3.21) 

The threshold level is also nicely coded in the LR tableaux: all elementary couplings 
have threshold level 1 and they all have a single column. We can then write directly that 

ko = ^columns = Ai + nn (3.22) 

and we recover the previous result. For an su(2) LR tableau, it is clear that the number of 
columns is given by this expression. More generally, for su(N), it is simple to check that 
the number of columns is simply 

N-l 

^columns = (A + \i + v,u>n-i) = Ai + nn (3.23) 

i=i 

where ujn-i is the iV — 1-th fundamental weight. 
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4. The generating function for su(S) fusion rules 



The su(3) tensor-product elementary couplings are: 



E l = 



E R 



Ep, 



E 7 



Es 



(4.1) 



Using the Kac-Walton formula, the threshold level of E\ is 1 and the corresponding 
fusion reads 

(4.2) 





E! : 


[0,1,0] x 


[0,0,1] D [1,0,0] 




Acting on E\ 


with (a n ,a m ;a n+m 


) n, m = 0, 1, 2 yields the elementary couplings 


Eq 


[1,0,0] x [1,0,0] 


D [1,0,0]: 


d 


(1,0,0,0,0,0,0,0) 


Ei 


[0,1,0] x [0,0,1] 


D [1,0,0]: 


dL 1 N 12 N 2 ^ 


(1,1,0,0,1,0,0,1) 


Ei 


[0,1,0] x [1,0,0] 


D [0,1,0]: 


dL l 


(1,1,0,0,0,0,0,0) 


Es 


[1,0,0] x [0,1,0] 


D [0,1,0]: 


dN n 


(1,0,0,1,0,0,0,0) 


E 4 


[0,0,1] x [0,1,0] 


D [1,0,0]: 


dL 2 N 13 


(1,0,1,0,0,1,0,0) 


E 5 


[0,0,1] x [1,0,0] 


D [0,0,1]: 


dL 2 


(1,0,1,0,0,0,0,0) 


Eq 


[1,0,0] x [0,0,1] 


D [0,0,1]: 


dN xl N 22 


(1,0,0,1,0,0,1,0) 


E 7 


[0,1,0] x [0,1,0] 


D [0,0,1]: 


dL x N 12 


(1,1,0,0,1,0,0,0) 


Eg 


[0,0,1] x [0,0,1] 


D [0,1,0]: 


dL 2 NuN 23 


(1,0,1,1,0,0,0,1) 


The last column is the vector e 


i with entries (k, X\, A2, nn, 


ni 2 ,nis,n 22 ,n 23 ). 



(4.3) 



By this 

procedure, we have thus recovered the affine extension of the 8 tensor-product elementary 
couplings and found an extra elementary coupling: Eq. 

To derive the fusion basis, we proceed as in the su(2) case. The set of variables here 



is 



(x , xi, • ■ ■ , x 7 ) = (k, Ai, A 2 , ran, 7112, rai 3 , n 2 2, n 23 ) 
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(4.4) 



and the matrix V (with columns written in the order Eq, • • • , Eg) reads 



V 



(4.5) 



/l 1111111 1^ 
1 1 1 
1 1 1 
1 1 1 
1 1 
000010000 
000000100 
Vo 1 l/ 

The reformulation of u T V > in terms of equalities by the introduction of appropriate 
nonnegative parameters reads: 

Uo = «o U0 + U2 = «5 

uo + u\ + U4 + U7 = a\ U0 + U3 + u 6 = a 6 

uo + u\ = «2 uo + u\ + 114 = a 7 (4.6) 

uo + u 3 = a s U0 + U2 + u 3 +u 7 = a 8 

Uo + U2 + U5 = «4 

We have 17 variables and 9 equations, hence 8 free variables. Let us choose them to be 
the ai except for a$. Solving for the dependent variables leads to 



Uo = 


a 


u 5 = 


— ao + ol\ + «3 + (X4 — ot-j 


Ul = 


— ao + a.2 


u 6 = 


-a 3 + a 6 


u 2 = 


—OL\ — Ct3 + OJ7 + OLg 


u 7 = 


Oil — CK7 


U3 = 


-a + a 3 


a 5 = 


CtQ — «i — «3 + «7 + ag 


114 = 


— a 2 + «7 







(4.7) 



The basis vectors of this system are obtained by setting one of the aj = 1 and all 
the others equal to (with the understanding the a 5 is excluded from this list of free 
variables). It appears more natural here to express them in their exponentiated version 
since a projection will be needed to extract the non-negative fundamental solutions. Denote 
by Ui the grading variable associated to Ui and by Ai those associated to 0^, the exponential 
form of the basis vectors reads 



So : U U^ 1 U 3 - 1 U^ 1 A A 5 
£1 : U^U^AxAt 
£2 : U 1 U^ 1 A 2 
£3 : U^UsUsU^AsAs 1 



£4 '■ U5A4 

£5 '■ UqAq 

£ G : U 2 U 4 U^ 1 U 7 - 1 A 5 A 7 
£ 7 : U 2 U^A 5 Ag 



(4.8) 
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To get the corresponding non- negative couplings, i.e., terms containing only non- negative 
powers of the Ai, we must keep only the non-negative powers of the £j. But this is not 
sufficient since negative powers of A5 can appear: we need to project the free generators 
of the non-negative S\ powers 

i=0 1 

to non-negative A$ powers, using, say the MacMahon algorithm (cf. the O projection 
described in section 3 of [1]). After the projection, all the variables Ai are set equal to 1. 
Here however, it is fairly easy to find out by inspection those non-negative combinations 
of the £i that have non-negative A5 terms. These are 

£0, £2, £4, £5, £6, £7 (4-10) 



together with 

£ £i : UoU^U^U^UrAoAi 
£ £ 3 : UqU^U^U^AqAs 
£l£ 6 : U4A1A7 



£i£ 7 : UtAtAa 

£ 3 £ 6 : UsUaU^U^AsAj (4.11) 
£ 3 £ 7 : U 3 U^A 3 A 8 



At this point, we set all Ai = 1. We have thus 12 elementary non-negative solutions 
and the corresponding inequalities are: 

Ai > n 12 A 2 > ni3 A 2 + n 12 > n 13 + n 23 

(4.12) 

nn > ri22 rin + nyi > ^22 + ^23 

and riij > (except for rin > which is implied by the others), which are the LR 
conditions for su(3). There are also three inequalities involving k: 

k - Ai - A 2 > n 22 

k — Ai - A 2 > nn - n 23 (4-13) 
k — Ai > riis + rin 

The set of inequalities (4.12) and (4.13) represents the su(3) fusion basis. 

Before we leave the analysis of the sit (3) case, let us return to the set of equations 
(4.7). The last equality gives a relation between different a*. Actually this relation signals 
a relation between different sums of columns of V. In other words, this signals a relation 
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between products of elementary couplings. Indeed, to link the last equality of (4.7) with 
such a relation, we recall that the labelling of the oti is that of the elementary couplings, 
which are the columns of V. Hence, the sought for relation is simply the product form of 
the equality with ctj i— > Ef. 



ai + Oi3 + = ao + 0L7 + as 1— > EiE 3 E 5 = EqE^Es 



(4.14) 



As there is only one relation, it is easy to find the generating function. Forbidding 

E-i_E 3 E 5 , we get [4] 



n c 1 -^) -1 



i=0 
i^l,3,5 



/ 



E-x 



(1-E 1 )(1-E 5 ) 
E3E5 



(4.15) 



+ 



(l-E 3 )(l-^i) (1-E 5 )(1-E 3 ) 



If instead, we decide to forbid EqE^Es, we would have 

/ 8 \ 



G' 



n a-^r 1 



i=0 



1 



(1-E )(1-E 7 ) 
E E 8 



(4.16) 



+ 



(l-E 7 )(l-E 8 ) (1-E 8 )(1-E ) 



and simple manipulations show that G\ = G' . An independent proof of this generating 
function is presented in Appendix A. 



Given the fusion basis, we can write down directly the threshold level to be 
ko = max(Ai + A 2 + nn - n 23 , Ai + A 2 + n 22 , Ai + nu + n 13 ) 



(4.17) 



This can also be extracted from the generating function as follows. A generic term of the 
su(3) generating function is (E = d) 



T? a T?b T? c I?d T? e T?f T?9 T?h 
&\ -^2-^3-^4 -^5-^6 -^7-^8 



with either a = 0, c = or e = 0. In all cases the threshold level is simply 



(4.18) 



ko = a + b + c + d + e + f + g + h 



(4.19) 
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In terms of the grading variables L; and Nij , the above generic term becomes 

dr+*oL$ +b+9 Li +e+h Nff f+h N£ 9 N? 3 N£ 2 N£ h 



(4.20) 



From this expression we read off the relation between the and the variables a, • • • , h. 
In each three cases (where one of a, c, e is zero), we can then solve for the sum ko = 
a + b + c + d + e + f + g + h. We find 



a = 0: k = Ai + A 2 + nn - 7i 2 3 
c = : k = Ai + A 2 + n 22 
e = 0: k = Xi + nu + n 13 



(4.21) 



This leads to the compact expression (4.17) for the sit (3) threshold level. This is easily 
checked to be equivalent to the formula given in [6,7] in terms of BZ triangle data (cf. 
section 7.1 of [1]): 



k = max{mi 3 + /zi + /u 2 , ™13 + ^i + ^2, Ji3 + Ai + A 2 } 
An explicit formula for the su(3) fusion coefficients is written down in [19]. 



(4.22) 



Notice that the threshold level is also simply encoded in the LR tableaux. Indeed, 
every elementary couplings has threshold level 1 and it corresponds to the number of 
columns except for Eg. This leads directly to the following formula for the threshold level 
of a general LR tableau 



k > k = # columns — #E$ = # columns — # 



(4.23) 



that is, ko is the number of columns minus the total number of Eg that we can take out 
of the tableau while preserving its LR character. Consider for instance: 









1 




1 


2 




1 















1 






1 




1 


2 










1 










_2 





(4.24) 



After the subtraction of one E 8 , the resulting tableau is not a LR tableau: counting from 
right to left, we find that a \2\ precedes the first |T|. Therefore, no E 8 can be removed 
and ko is given by the number of columns which is 4. 
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5. The sp(A) generating function 

We first recall some results obtained in [1] . The appropriate basis for the description 
of sp(4) tensor products reads [20]: 

Ai > p iii > q 

A2 > ri/2 hi > q + v\ - r 2 

\ 2 >r 1 /2 + q-p Hi>P + r 1 -r 2 (5.1) 

A 2 >r 2 /2 + g-p H2>r 2 /2 

v\ = r 2 - T\ - 2p + Ai + Hi h> 2 = p - q - r 2 + A 2 + /U 2 

together with p, q e N and r» G 2N for i = 1, 2. A proper set of variables for a complete 
description of a particular tensor-product coupling is thus 

{\ 1 ,\ 2 ,Hi,V2,r 1 ,r 2 ,p,q} (5.2) 

(notice the absence of the Vi Dynkin labels). Let the corresponding grading variables be 

{L u L 2 , M u M 2 , i? 2 , P, Q} (5.3) 

The list of elementary coupling with their grading description is: 



An 


(0,0) ® 


(1,0) 


D 


(1,0) 


M x 


A 2 : 


(1,0) ® 


(0,0) 


D 


(1,0) 


£1 


A 3 : 


(1,0) <g> 


(1,0) 


D 


(0,0) 


L X M X PQ 


B l : 


(0,0) ® 


(0,1) 


D 


(0,1) 


M 2 


B 2 : 


(0,1)® 


(0,0) 


D 


(0,1) 


L 2 


B 3 : 


(0,1)® 


(0,1) 


D 


(0,0) 


l 2 m 2 r\r 


Ci : 


(0,1)® 


(1,0) 


D 


(1,0) 


L 2 MiQ 


C 2 : 


(1,0) 


(0,1) 


D 


(1,0) 


L X M 2 R\P 


<^ 3 : 


(1,0) 


(1,0) 


D 


(0,1) 




£>l : 


(2,0) 


(0,1) 


D 


(0,1) 


l\m 2 r\p 


£>2 : 


(0,1)® 


(2,0) 


D 


(0,1) 


L 2 MlR\ 


D 3 : 


(0,1)® 


(0,1) 


D 


(2,0) 


L 2 M 2 R\ 



(5.4) 
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(5.5) 



The relation between elementary couplings are generated by 

C 1 C 2 =A 3 D 3 , C 2 C 3 = A l D 1 C 3 C 1 = A 1 A 3 B 2 

D 1 D 2 = B 3 Ci D 2 D 3 = A\B 2 B 3 D 1 D 3 = B 2 C\ 

C\D\ = A 3 B 2 C 2 C 2 D 2 = A\B 3 C 3 C 3 D 3 = A\B 2 C 2 

To find the fusion elementary couplings, we start by computing the threshold level 
of A\ by the Kac- Walton formula. It is found to be 1. The corresponding level- 1 fusion, 
denoted A 1: is thus 

[1,0,0] x [0,1,0] D [0,1,0] (5.6) 

We can act on it with the four pairs 

(A A') = {(1, 1), (a, a), (a, 1), (1, a)} (5.7) 

We obtain in this way two copies of A\ and two copies of C\, the level-1 extension of C\. 
Similarly, A 2 , and A 3 are found to have level 1 and this implies the same result for C 2 , C 3 . 
Bi is also found to have threshold level 1. Acting on it with the above sequence of outer 
automorphisms leads successively to B\, B 2 , B 3 and a new coupling, E : 

E : [1,0,0] x [1,0,0] D [1,0,0] (5.8) 

Finally, Di,D 2 and D 3 have threshold level 2 and they are all fixed with respect to the 
action of the outer-automorphism group. The set {Ai, Bi, Ci, Di, Eq] is thus our candidate 
complete set of fusion elementary couplings, whose explicit expression in terms of grading 
variables is read from their tensor-product relative with the addition of an appropriate 
factors of d. 



Having obtained the fusion elementary couplings, we now work out the corresponding 
fusion basis. Introduce the set of variables 



x 8 ) 



(k,X 1 ,X 2 ,ii 1 ,ii 2 ,r 1 ,r 2 ,p, q) 



(5.9) 



This fixes the ordering of the rows of V. The matrix V is built from the columns which 
form the different elementary couplings in the order Eq, Ai, ■ ■ ■ , D 3 : 



(1 11111111122 2\ 
0011000011200 
000001110001 

10 10 2 



V 



1 

0101000101020 
0000101010101 
0000002000020 
2 







2 2 



00010000112 



2 




(5.10) 



\0 00100010000 0/ 
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The transcription of the inequalities u T V > into the equalities u T V = a T takes the 
following form 

Uq = « Uq + U2+U3 + u 8 = «7 

"^0 + ^3 = OL\ Uq + U\ + U 4 + 2tt6 + ^7 = «8 

«0 + Wl = «2 «o + Ui + Us + U 7 = Ctg 

u + ui + u 3 + u 7 + u 8 = a 3 2u + 2wi + u 4 + 2u 6 + 2u 7 = a w (5.11) 

uq + u 4 = a 4 2u + u 2 + 2u 3 + 2u$ = an 

u + u 2 = a 5 2uq + u 2 + u 4 + 2u 6 = ai 2 

uo + u 2 + u 4 + 2-U5 + 2uq = a& 
Solving for the dependent variables Ui, aj, i = 0, • • • , 8 and j = 6, 7, 8, 9 gives 

uo = cto U7 = \{olq — 2a 2 + «5 + aio — ai 2 ) 

ui = —a + a 2 u 8 = \{cxq — 2ai + 2a 3 — a 5 — aio + 012) 

u 2 = — «o + «5 «6 = — 2«i — a 5 + an + «i2 

us = -a + ai a 7 = \{-a + 2a 3 + a 5 - a w + ai 2 ) (5.12) 

u 4 = — a + a 4 a 8 = |(— cto — «5 + «io + 0=12) 

w.5 = ^(«o — 2cki — «5 + «n) a 9 = 7j(— a + 2cti + a 5 + a w — ai 2 ) 

u 6 = |(-Q!4 - «5 + "12) 
As usual, the basis vectors of this system are obtained by setting one of the ai = \ and 
all the others equal to 0, excluding a$, ■ ■ ■ , ag. We will give their exponentiated version, 
where as before, we denote by Ui the grading variable associated to ui and by A% those 
associated to af. 

51 : UsU^U^AiAfAg 

5 2 : UiU 7 ~ 1 A 2 
Ss : U 8 A 3 A 7 

5 4 : U 4 U~ 1/2 A 4 (5.13) 

5 5 : t/^-^Ve- 1 /^, 1 / 2 ^- 1 / 2 ^- 1 ^ 72 ^ 172 ^ 72 
<?6 : U 7 ^ 2 U 8 Aj 1 ^ 2 A^ 2 A^ 2 Aio 

S 7 : £/ 5 1/2 A-4ii 

£ 8 : t/ 6 1/2 t/ 7 - l72 ^ 1/2 A^ 7 72 ^ 72 A- l72 ^2 
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Next we keep only those combinations of the 8i that contain only non-negative integer 
powers of the A4. This projection is not so simple to work out by inspection. We thus 
need to use a more systematic procedure: 

Consider a general expansion of the form ]^[^ (1 — and in a generic term of the form 
[li^S l e t us collect the number of Ai factors (denote by a; their exponents). Of course 
we are only interested in those Ai that appear with negative powers, namely i = 6, 7, 8, 9. 
Their powers can be read off from the Ai in (5.13) and this yields the following expressions: 

a 6 = —2e[ - €5 + e' 7 + eg > 

2a 7 = -e' + 2e' 3 + e' 5 -e' 6 + e' 8 >0 

(5.14) 

2a 8 = -e' -e' 5 + e' 6 + e' 8 >0 

2a 9 = -e' Q + 2e[ + e' 5 + 4 - e 8 > 

(These equations should be compared with the last four of (5.12), with ai — > e[, i < 5 
and ai — > e-_ 4 for i > 10). We then look for the elementary solutions of this system of 
inequalities. There are 4 elementary solutions with 7^ 0. Their grading reformulation 
reads 

SqSiSqEI 80818788 8o8s8q8s 8q8^8q (5.15) 

Denote their vector-reformulation respectively as with i = 0, 1, 2, 3, then the conditions 
Gi x > yield, in the above order 

k > Ai + A 2 + [12 + ri/2 - r 2 
k>X 1 + X 2 + iJ2- r 2 /2 
k > Ai + fii + fi 2 - V 
k > Ai + A 2 + fii + \i 2 - V - 1 - r i/2 

The other elementary solutions are 

82, 8s, 84, 87, 858s, 8q8%, 8\8\, 8187, 
£38%, 858587, 858^88, 8\8q8%, 8x8^878% 

and the resulting inequalities reproduce the whole set of BZ inequalities (5.1) with the 
positivity requirement on ri,p and q (together with /ji > q + \ (ri — r 2 ) which is implied 
by the other ones). 
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(5.16) 



(5.17) 



Let us return to the last four equations in (5.12). As mentioned in connection to the 
sit (3) case, they indicate the 'basic relations': the correspondence between the aij and the 
elementary couplings being fixed by the ordering of the columns of V (e.g., a 3 h- > As and 
«7 | — C\). The relations correspond then respectively to 

A\B 2 Bs = D 2 D 3 E dfD 1 = A 2 3 B 2 D 3 

(5.18) 

E B 2 Cl = D x Ds E Q ClDs = A 2 1 B 2 D 1 

The first and third relations appear in the list (5.5). All other linear relations in the set 
(5.5) can be obtained from products of the above four, allowing for the cancellations of 
common factors. For instance, consider the product of the left factors of the second and 
third relations; equating this with the product of the right factors yields 

E%B 2 d 2 1 C 2 2 D 1 = A\B 2 DiL>1 (5.19) 

Cancelling the B 2 D\ terms and taking the square root gives 

%A C 2 = AsDs (5.20) 

which is the affine extension of the relation C\C 2 = A3.D3. All other linear relations can 
be obtained in a similar way. 

We can write the sp(4) generating function in the compact form 



G = E B 1 B 2 Bs[A 1 A 2 A 3 C 1 C 2 Cs(l - A x AsB 2 ) + D 1 D 1 A 2 A 3 C 2 C 3 
+ D 3 D 3 A 1 A 2 C 1 C 2 + D^MMAsC^^l - A^A)] 



(5.21) 



where Q is defined as 

Q=— (5.22) 
1-Q 

This can be re-expressed under a manifestly positive form as follows 

G = EqA\A 2 B\B 2 B 3 C\C 2 C 3 + EqA 3 A 2 A 3 BiB 2 B 3 CiC 2 C 3 
+ EoAtAsA^AsB^sC^Cs + E^MMB^BsC^^ 
+ EqDsA\A 2 BiB 2 BsC ±C 2 Ds + EqD 2 A\A 2 BiB 2 BsC \CsD 2 
+ EoA 3 D 2 A 2 A 3 B 1 B 2 B 3 C 1 C 3 D 2 + E q A x A z D 2 A~Jl 2 A~ 3 B{B z C{C z D 2 
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(5.23) 



We should stress that this is essentially a new result. A generating function for sp(4) fusion 
rules was given in [21]; the approach, however, was ad hoc and the result was not related 
to any known basis. 

As before the information concerning the threshold level that can be deduced from the 
fusion basis inequalities (5.16) can also be obtained directly from the generating function. 
A generic term of the sp(A) generating function (5.23) reads 

d a A1A b 2 A c 3 BfB^B(Cfd^C l 3 D{D^D l 3 (5.24) 

Its threshold level is (since all these factors have a single power of d except for the three 
Di = d 2 Di): 

k = a + b + c + d + e + f + g + h + i + 2j + 2k + 2l (5.25) 

Now express the elementary couplings in terms of dummy variables {Li, L 2 , Mi, M 2 , R\, R2, P, Q} 
whose exponent are the BZ basis data, respectively {Ai, A 2 , /zi, H2, r i, r 2,P, q}'- 



A 1 = Mi A 2 = L 1 A 3 = L 1 M 1 PQ 

Bi = M 2 B 2 = L2 B3 = L2M2R1R2 

d = L 2 M 1 Q C 2 = L X M2R\P C 3 = L X M X P 

D 1 = LlM 2 RlP 2 D 2 = L 2 MlR\ D 3 = L 2 M 2 Rj 



(5.26) 



Next, consider each term of the generating function (5.23) and solve for k in terms of the 
basis variables. Surprisingly there are only four different formulas for k . The expressions 
corresponding to the different terms of (5.23) are: 

terms 1, 5, 6 : ko = Ai + A 2 + ^1 + ^ 2 - p - q - n/2 
terms 3,8: ko = Ai + [i\ + /U 2 — p 

(5.27) 

term 7: k = X 1 + A 2 + /i 2 + r 1 /2 - r 2 

terms 2, 4 : k = Xi + A 2 + /j 2 - r 2 /2 

Therefore, the threshold formula is the maximum value of these four values or equivalently 

fc = Ai + A 2 + A*i + A*2 - min(p + q + n/2, A 2 + p, fii - n/2 + r 2 , \i\ + r 2 /2) (5.28) 

Notice that by rewriting k > ko, we recover from (5.27) the 4 inequalities (5.16). 

26 



The system of inequalities (5.1) can be transformed into a system of equations by 
setting n/2 = si and T2/2 = s 2 and introducing the integers ai (cf. section 7.5 of [1]): 



Ai 


= V + «i 


V 2 = <X4 + as 


A 2 


= si + a 2 


a 2 +p = a 3 +q 


m 


= q + a 5 


CI3 + Si = CI4 + S2 




= s 2 + a 8 


a 5 + 2s 2 = a 6 + 2si 


V\ 


= cii + a-? 


a 6 + q = a 7 + p 



As shown in [1], this leads to the following diamond-type graphical representation of the 
tensor product: 




Dotted lines relate those two points that compose the label indicated beside it and 
opposite continuous lines are constrained to be equal, with the length of a line being defined 
as the sum of its extremal points except for the lines delimited by the points (ciq, si) and 
(05,^2) where the point Si is counted twice (the little bar besides si and s 2 being a 
reminder of this particularity). For those lines, the constraint reads a 6 + 2si = a 5 + 2s 2 . 
Given a triple sp(4) product, the number of such diamonds that can be drawn with only 
non-negative entries gives its multiplicity. 
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In terms of these data, the expression for the threshold level (5.27) look somewhat 
more symmetrical: the four expressions in (5.27) correspond respectively to the following 
terms: 



k = a 1 +a 8 + max{a 4 + a 7 + si, a 5 + q + s 2 , a 4 + q + si, a 4 + q + s 2 } (5.30) 



6. The sit (4) generating function 

Written directly in terms of LR tableaux, the sit(4) elementary solutions are: 



At = 



a 3 = u, b x = 



B 2 = 



Ci = a, c 2 



c 3 



, ^2 



#3 = 



> ^3 



(6.1) 



and 







1 


D 1 = 


2 






1 









1 






1 


£> 2 - 




2 


, ^3" 








_3 






_2_ 





En 





1 




2 


1 




1 





(6.2) 



The relations are [22,23]: 



D,D k = CiEi 



DjDk = BiCjC k EiEj = B k D k D k 



DiEi = CjB k D k = BjDjC k 

with i, j, k a cyclic permutation of 1, 2, 3. 



(6.3) 



Consider now the construction of the set of fusion elementary couplings using outer- 
automorphism completeness. Start with A\ : (0, 0, 0) ® (0, 0, 1) D (0, 0, 1), this has thresh- 
old level 1. Acting on it with 



where 



(A, A') = (a n , a m ), n,m = 0,1, 2, 3 

o[Ao, Ai A 2 , A3] = [A3, Ao Ai, A 2 ] 
28 



(6.4) 



(6.5) 



we generate the affine extension of the whole set A^B^C^D^^D'^ which thus all have 
threshold level 1, together with the scalar coupling E = [1, 0, 0, 0] x [1, 0, 0, 0] D [1, 0, 0, 0]. 
Finally, the affine extension of Ei arises first at level 2: [0, 1, 0, 1] x [1, 0, 1, 0] D [1, 0, 1, 0]. 
The three weights in this coupling are fixed under the action of A = a 2 . Hence, we need 
only to consider 

(A A') = {(1, 1), (a, a), (a, 1), (1, a)} (6.6) 

and this leads respectively to Ei, E2, E%, which all have ko = 2 and a new elementary 
coupling F = [0, 1, 0, 1] x [0, 1, 0, 1] D [0, 1, 0, 1] (first discovered in [7]). Notice that at the 
level of tensor-products, F is a composite product C\C 2 C 3 . But if it were still composite 
for fusions, it would necessarily have level 3 since k (CiC 2 Cs) = 3. This is the reason why 
F must be regarded as a new elementary coupling. 



The whole set of fusion elementary couplings is: 



M = 


[1,0,0,0] x [0,0,0,1] 


D 


[0,0,0,1] 


D[ = 


[0,0,1,0] 


x 


[0,1,0,0] 


D 


[0,0,0,1] 


A 2 = 


[0,0,0,1] x [0,1,0,0] 


D 


[1,0,0,0] 


D' 2 = 


[0,1,0,0] 


X 


[0,1,0,0] 


D 


[0,0,1,0] 


M = 


[0,1,0,0] x [1,0,0,0] 


D 


[0,1,0,0] 


D' 3 = 


[0,1,0,0] 


X 


[0,0,1,0] 


D 


[0,0,0,1] 


Bi = 


[0,0,0,0] x [0,0,1,0] 


D 


[0,0,1,0] 


Di = 


[0,0,1,0] 


X 


[0,0,0,1] 


D 


[0,1,0,0] 


B 2 = 


[0,0,1,0] x [0,0,1,0] 


D 


[1,0,0,0] 


D 2 = 


[0,0,0,1] 


X 


[0,0,0,1] 


D 


[0,0,1,0] 


B 3 = 


[0,0,1,0] x [1,0,0,0] 


D 


[0,0,1,0] 


D 3 = 


[0,0,0,1] 


X 


[0,0,1,0] 


D 


[0,1,0,0] 


C 1 = 


[1,0,0,0] x [0,1,0,0] 


D 


[0,1,0,0] 


Ei = 


[0,1,0,1] 


X 


[1,0,1,0] 


D 


[1,0,1,0] 


c 2 = 


[0,1,0,0] x [0,0,0,1] 


D 


[1,0,0,0] 


E 2 = 


[1,0,1,0] 


X 


[1,0,1,0] 


D 


[0,1,0,1] 


C 3 = 


[0,0,0,1] x [1,0,0,0] 


D 


[0,0,0,1] 


E3 = 


[1,0,1,0] 


X 


[0,1,0,1] 


D 


[1,0,1,0] 



together with two couplings that have no elementary finite relative: 



£0 = [1, 0, 0, 0] x [1, 0, 0, 0] D [1, 0, 0, 0] F = [0, 1, 0, 1] x [0, 1, 0, 1] D [0, 1, 0, 1] (6.8) 

The tensor-product relations are modified by the appropriate insertions of d or Eq 
factors in order to put them at the same threshold level: 

E D-Dk = CiEi E DjD k = BiCjCk EiE~ = E BkDkD k 

3 _ ; ,J ^ ^ ^ ^ ^ (6.9) 

DiEi = CjBkDk = BjDjCk EqF = C\C 2 C 3 

with i, j, k a cyclic permutation of 1, 2, 3. 
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To get the su(4) basis, we first write down the V matrix, whose columns are the vec- 
torial transcription of the fusion elementary couplings written in terms of the grading vari- 
ables. The column ordering corresponds to -Bo, Ai, Bi, Ci, D[, Di, Ei, F with i = 1,2, 3. The 
rows are labelled by the LR variables (k, \\, A2, A3, ri\\,n\2, ri\z, ri\±, 7122, ^23, ^24, ^33, ^34)- 
The matrix V is thus 



V 



Z 1 



1 1 













111111 

1 1 

1 1 



111111 

110 







1 1 



1 







100000010000 



1 1 

100100100000111 
000010011000 








2 

1 



1 



1 



2 2 2 

1 

1 1 



1 1 



010000100000001 



0010000000000000000 



010010000000 
000000001000 





1 1 



1 









100000000011 





1 



1 







0100000000000000000 





1 
1 



1 







1 






\ 



\o 0000000100001100011 / 



(6.10) 



With u = (uo, • • • , U12), the equations u T V > can be transformed into equalities by 
introducing the variables af. 



Uq = a 

U + U 4 + U 8 + Un = OL\ 

uq + u 3 + u 7 = a 2 

Uq + U\ = tt 3 

Uq + U4 + Us = OL4 

UQ + U 2 + Uq + Uio = a 5 

Uq + U2 = OLQ 

Uq + U4 = Ct7 

Uq + U\ + U 5 + U 9 + U\2 = «8 

u +u 3 = a 9 



U0+U2+U6 = Ctio 

u + ux + u 5 = an 

U + Ui + U 5 + U 9 = «12 

U + U 2 + U4 + U 9 + U12 = «13 

U + U 3 + U4 + U 8 + U\2 = CK14 

U0+U3 + U4 + Uio = «15 

2-u + Ui + U3 + U 5 + Uio = "16 

2u + u 2 + U4 + u 9 = an 

2u + U 2 + U4 + U 6 + U S + U12 = «18 

2u + Ui + U 3 + U4 + U 5 + U 9 + U12 = aig 



(6.11) 



We have 13 free variables; let us choose them to be the ai for % = 0, • • • , 12. Solving 
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for the dependent variables leads to 



u 


= «0 


W 5 = 


-a 3 + an 


MlO = «5 - "10 


1 


= — Ctn + 


Uq = 


— CtR + am 


tin = Cti — Ct4 


U 2 


= — a + «6 


U 7 = 


1X2 tx 9 


u 12 — "8 "12 




= — «o + «9 


U 8 = 






U4 


= — Cto + «7 


Ug = 


— CKH + «12 




with 










«13 = 


— a + «6 + «7 + «8 — 


an «i 7 


= a 6 + a 7 — «n + «12 


"14 = 


— «0 + OC4 + Ct8 + Ct9 — 


«12 "18 


= a 4 + a 8 + «io — «i2 


«15 = 


— ttQ + «5 + OL 7 + Ctg — 


ttlO «19 


= — tto + «7 + «8 + «9 



(6.12) 



(6.13) 

— CK10 «19 = — CtO + «7 + «8 + «9 

«16 = «5 + "9 — «10 + "11 

Now, by setting successively = 1 for i = 0, • • • , 12 and the others equal to 0, we 
generate the following set of basis vectors: 

So = dN^L^L^L^AoA^A^A^A^ 1 £ 6 = AT" 1 L 2 A 6 A 13 A 17 

£ l = N 33 At £ 7 = N 11 N^L 2 A 7 A 13 Ai b A 17 A 1Q 

£ 2 = N 1A A 2 £ 8 = N 34 A 8 A 13 A 14 A 18 A 19 

£3 = N^LiAj, £9 = -^I4 1 -^3^9^14-4l5^16^4l9 

£ 4 = N^-N 33 A A A 1A A 18 £ w = N 13 N^ A w A^A^A 18 

£5 = N^AsAisAw £11 = A^i2A^ r 2 3 1 ^4ii^4 13 1 ^li 6 ^l 17 1 

£12 = N 23 N^A 12 A^A 17 A^ 

(6.14) 

We must now look for those combinations that contain only non-negative powers of 
the Ai. Since each £j contains at least one positive power of A4, these must be obtained 
from positive combinations of the £j. To find them, it is convenient to proceed as in the 
analysis of sp(4) . Denote by ai the number of Ai factors in a general term Yl £^ of the 
free expansion of the £{ in non-negative powers we get ( equivalent ly, we can read off the 
Ai from (6.14)): 

Ol3 = -Co + e 6 + e 7 + e 8 ~ e n a 17 = 4 + 4 - e il + e 12 

a i4 — —^0 + e 4 + e 8 -+- e 9 — e 12 Ois — e 4 + e 8 + e io ~~ e i2 
ais = — e' + e' 5 + e' 7 + e' Q — e' 10 aig = — e + e 7 + e§ + e 9 

ai6 = 4 + 4 _ e io + e n 
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(6.15) 



These relations are to be compared with (6.13). We thus look for elementary solutions of 
the system > for e\ non-negative integers. 

The full list of composites that involve 8q - these are those that generate /c-dependent 
constraints - is 

808487, 808389, 808589, 808983, £08789, 80858s, 

2 (6-16) 

808783811, 808789812, 808789810, 80858381284 80878389 

The constraints are (£* specifies the vector ej such that the inequality is ej x > 0) : 



808487 : 


k > Ai 


+ A 2 + A 3 + n 33 


808389 : 


k > Ai 


+ A 2 + + ^14 - n 3 4 


808589 : 


k > Ai 


+ nn + n i3 + ni4 


808783 : 


k > Ai 


+ A 2 + A 3 + n 22 — n 34 


808789 : 


> Ai 


+ A 2 + 77-14 + n 22 


808583 : 


> Ai 


+ A 2 + A 3 + nn - n 24 - n 34 


808783811 : 


> Ai 


+ A 2 + A 3 - ni2 + n 22 + n 23 - n 34 


808789810 : 


> Ai 


+ A 2 + ni 4 - ni3 + n 22 + n 24 


808389812 : 


k > Ai 


+ A 2 + nn + n i4 - n 23 


80878389 : 


2/c > 2Ai + 2A 2 + A 3 + n i4 + n 2 2 + nn - n 34 



(6.17) 



When these inequalities are re-expressed in terms of BZ triangle data, they reproduce the 
threshold formula presented in [7]. 

The £cr m dependent elementary solutions, namely 

81, 82, 83, 84, 85, 85, 87, 83, 89 

85811, 84812, 87811, 83812, 89810, 85810, (6.18) 

^8^11^12, 89810812, ^7^10^11 

yield the standard LR inequalities: 

Ai > n i2 nn > n 2 2 

A 2 > n i3 nn +n i2 > n 22 + n 23 

A2 + n i2 > ni3 + n 2 3 nn + n i2 + n i3 > n 2 2 + n 2 3 + n 2 4 

(6.19) 

A3 > n M n 2 2 > n 33 

A3 + n i3 > ni4 + n 2 4 n 2 2 + n 2 3 > n 33 + n 34 

A3 + n i3 + n 23 > n M + n 24 + n 34 
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and riij > 0, except for ri\\ > and n 2 2 > which are implied by the above equations. 

As in the sp(4) case, we can check that the relations (6.13) code the 'basic linear 
relations' of the model. Indeed, the 7 relations read from (6.13) are 

EqDiD' 2 = B 3 C\C 2 EiD[ = B 2 C 3 D' 2 

J J J J J J J J ( 6 - 20 ) 

EqD 3 D[ = B2C1C3 E 3 D' 3 = BiC 3 D[ 
EqF = C1C2C3 

and these are the generators of all the su(4) linear relations. 

In order to construct the sit (4) generating function, we must choose a term ordering. 
We fix the ordering as follows: 

{Li, L 2 , L 3 , N lh N 12 , N 13 , N 14 , iV 22 , N 23 , N 24 , N 33 , iV 34 , d, 

(6.21) 

E±, E 2 , E 3 , Bi, B 2 , B 3l Ci, C 2 , C 3 , Ai, A 2 , A 3 , Di, D 2 , D 3 , D±, D 2 , D 3 , E , F} 

Grobner basis methods yield the forbidden products: 

{EiEj, D i E i , DiEi, CiEi, BiCjCk, B 2 C 3 D 1 D 2 , B 1 C 3 D 1 D 2: B 1 C 2 D 3 D 1 , 

_ , 1_ . (6.22) 

E B 1 D 1 D 3 D 1 D 2 , FEi,FBi, C X C 2 C 3 , FC X C 2 C 3 } 

The different terms of the generating function are fully specified by their denominator (the 
numerators are introduced to avoid over-counting): 

EoMA^ADAd'^D'^d^, E A 1 A 2 A 3 B 2 B 3 d 2 d 3 D 2 D 3 D' 2 3' 3 E 1 F, 

E AiA 2 A 3 BiB 2 B 3 C 2 DiD 2 d' 1 d' 2 E 3 F, EoA^AsB^BsC^Ad'^EsF, 

EoAiAiAsBiC&DADsD'^D's, %A 1 A 2 A z B 1 %C 1 C 2 D 1 D 2 D[D 2 D^ 

E A 1 A 2 A 3 B 1 B2C 1 d 2 D 1 D 2 D 3 D 2 D 3 , EoA 1 A 2 A 3 B 3 C 1 C 3 D 1 D 2 D 3 d[d 2 D 3 , 

E A 1 A 2 A 3 BiB 3 diC 3 D 2 D 3 D 1 D 2 D' 3 , EoA 1 A 2 A 3 B 1 B 3 CiC 3 D 1 D 2 D 1 D 2 d' 3 , 

(6.23) 

hA l A 2 A 3 B l %C 1 C 3 D l D 2 %D 2 D 3 , ^MA^B^C^dA^D^ 

A 1 A 2 A 3 B 1 B 3 C 1 d 3 D 1 D 2 D 3 D 1 D 2 D 3 , A 1 A 2 A 3 B 1 B 2 B 3 C 1 D 1 D 2 D 3 D 1 D 2 D 
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E A 1 A 2 A 3 B 2 B 3 C 1 D 1 D 2 D 3 D 1 d' 2 D 3 , EoA 1 A 2 A 3 B 1 B 2 B 3 CiD 2 D 3 D 1 D 2 D 3 , 

e a 1 a 2 a 3 b 1 b 2 b 3 c 1 d 1 d 2 d' 1 d' 2 d' 3i e a 1 a 2 a 3 b 1 b 2 b 3 c 1 d 1 d 2 d 3 d' 2 d 3 , 

EoA^A^B^^dADsD^ 
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Finally, note that the expression of threshold levels in terms of tableau data, that is, 
the analogue of the su{2, 3) formulas written previously is clearly 

k = #columns - max{#D 1 + #D 2 + #D 3 + #CiC 2 C 3 } (6.24) 

since the D^s have level 1 but two columns and C1C2C3 has level 2 and three columns 
(corresponding to the F fusion coupling.) 

7. Conclusion and open problems 

We have obtained the fusion generating function for su(3, 4) and sp(4) using the con- 
jectural existence of a fusion basis. In the su(3) case a first-principle derivation (presented 
in Appendix A) provides an independent proof of the results, thus a partial confirmation 
of the conjectures and the correctness of the underlying fusion basis. Moreover, different 
tests of the sp(4) and sit (4) generating functions, presented in Appendix A, also support 
our conjectures and the fusion basis constructions. En passant, we point out that the 
search for the complete su(N) level-rank symmetric function introduced in Appendix A is 
a quest that deserves further studies. 

Although the theme of this paper is the construction of fusion generating functions, 
our most important result is the unravelling of the fusion basis concept, for which we 
have provided concrete examples. The main open problem is to find a fundamental and 
Lie algebraic way of deriving the fusion basis (analogous to the Berenstein-Zelevinsky 
conjectures [20]). We observe that the number of fc-type inequalities increases rather 
quickly with the rank of the algebra: 1 for stt(2), 3 for stt(3), 4 for sp(4) and 10 for su(4). 
More specifically we would like to find arguments to justify the homogeneity property (on 
the other hand, the linearity appears to be a generic property, a direct consequence of the 
Kac- Walton algorithm). 

With regard to the automorphism completeness conjecture we note that for simplicity 
(and because the discussion is to a large extent devoted to su(N) for which the outer- 
automorphism group is rather large) we have focused on the outer-automorphism group 
as the essential symmetry. It is natural to extend the conjecture to the full symmetry 
group of fusion coefficients. However, we should stress is that the outer-automorphism 
conjecture is just a convenient tool. If the conjecture (or its natural extension to the full 
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fusion symmetry group) turns out to be wrong, there are other avenues that could yield 
the complete set of fusion elementary couplings. 

In the present work, the only information on fusion data that has been extracted, out 
of the fusion basis or the fusion generating function, is the expression for the threshold level 
in terms of the basis variables. But there are certainly more data that can be lifted. For 
example, given a triple product with multiplicity m, to which there correspond m values of 
the threshold levels, we could ask for the expression, in terms of the Dynkin labels, of the 
minimum and maximum values of kg. It is easy to write down some explicit expressions 
for particular fusion coefficients. 

The reformulation of the problem of computing fusion rules in terms of a fusion basis 
solves, in principle, the quest for a combinatorial method since it reduces a fusion compu- 
tation to solving inequalities. But we expect that we have not found an optimal solution 
to the quest for an efficient combinatorial description. 
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Appendix A. Independent verifications of the fusion generating functions 

The su(3) fusion generating function is not presented here for the first time; it ap- 
peared originally in [4]. A sketch of its proof was presented in [24] without details. In 
this section we present a complete proof of the su(3) generating function for fusion rules; 
in addition, we describe some independent checks confirming the validity of the su(4) and 
sp(4) fusion generating functions given in sections 6.3 and 7.3. The first check that we 
present uses Giambelli-type formulas. These can be viewed as equalities of corresponding 
expressions in the character rings. Since the su(n) and sp(n) fusion rings are quotients 
of the classical character rings (see [25] and references therein), these formulas continue 
to hold for fusion products. For stt(4), we present another non-trivial check based on a 
level-rank duality argument. 
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A.l. Determinantal formula and the 'composition' method: deriving the su(3) generating 
function for tensor products 



The Giambelli formula, or more generally, determinantal formulae which give expres- 
sions for group characters as determinants, provide another method for calculating fusion 
generating functions in terms of simpler generating functions. This uses the technique of 
'composition' of generating functions described previously in section 2.3 of [1]. 

The sit (3) Giambelli formula expresses a general representation in terms of a difference 
of products of representations with a single non-zero Dynkin label, i.e., 

(Ai, A 2 ) = (Ai + A 2 , 0) <g> (A 2 , 0) - (Ai + A 2 + 1, 0) <g> (A 2 - 1, 0) (A.l) 

This can be rewritten in determinantal form as follows 

(\ A )-d P ff ( A i + A 2>°) (A 2 -1,0)\ 

(Al ' A2) - det Ui + A 2 + l,0) (A 2 ,0) ) (A ' 2) 

Consider first the generating function Gi(Li, L 2 , M 1? i? l5 R 2 ) which is the generating 
function for products of the form: (Ai, A 2 ) ® (/xi, 0). Its explicit form is 

Gl = (1 - LiiVi)(l - L 2 N 2 )(1 - L 2 M 1 )(1 - M 1 N 1 )(1 - L 1 M 1 N 2 ) (A ' 3) 

It is obtained by setting M 2 = in the complete tensor-product generating function 
(cf. section 2.5 in [1]). Our point here is not to re-derive G± from first principles but 
simply to show how we can reconstruct the complete generating function out of the partial 
information contained in G\. In the fusion case, we will indicate how the analogue of G\ 
can be obtained, preventing the argument from being circular. 

From two copies of G\ we form the composite generating function G 2 : 

G 2 (Li, L 2 , Mi, M 2 , iVx, N 2 ) = Q G^L^ M 1? R t , R 2 )G 1 (R^\ R 2 \ M 2 , N u N 2 ) (A.4) 

which is the generating function for products of the form 

(Ai,A 2 ) ® (//i,0) ® Q/2,0) (A.5) 

Note that the generating function for products 

(Ai,A 2 )(g) (ah + 1,0)® (// 2 -l,0) (A.6) 
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is M 2 M 1 G 2 and so, by (A.l), the generating function for products (Ai, A 2 ) <8> a* 2 ) is: 

Mi , 

G 3 = Q (G 2 - M 2 M~ 1 G 2 ) (A.7) 

The coefficient of M^M^ 2 is the multiplicity of the representation with Dynkin labels 
(Hi — [i 2 , n 2 ) in the product 

(Ax, A 2 ) <g> [(//!, 0) ® (/x 2 , 0) - + 1,0)® (/x 2 - 1, 0)] (A.8) 

To change to variables which carry the Dynkin labels we make the substitution M 2 1— > 
M 2 M 1 _1 , so that Mi now carries the first Dynkin label. This introduces negative powers 
of Mi, corresponding to products (A.8) with ^1 < /U 2 , which are not required. So we must 
keep only non-negative degree terms in Mi to obtain the final generating function. Denote 
the resulting expression as G^{L\,L 2l Mi, M 2 , Ni,N 2 ); it reads 

(1 - LiL 2 MiM 2 NiN 2 ) 

G4 = 



(1 - LiNi)(l - LxM 2 )(l - L 2 Mi)(l - L 2 N 2 ) ^ g) 



x 



(1 - M 2 N 2 )(1 - MiNi)(l - LiMiN 2 )(l - L 2 M 2 Ni) 
which is the usual form of the su(3) generating function (cf. section 2.5 of [1]). 



A. 2. Extension of the determinantal formula methods to fusion rules: the su(3) case 

The starting point for the derivation of the su(3) fusion generating function is the 
generating function for fusions of the form 

[k - Ai - A 2 , Ai, A 2 ] x [k - /xi,//i, 0] (A.10) 

These fusions are known explicitly and the information on their fusion coefficients can be 
lifted to the following generating function [26] 

Fi(d, Li, L 2 , Mi, Ni, N 2 ) = 

1 (A.11) 

(1 - d)(l - dLiNi)(l - dL 2 N 2 ){\ - dL 2 Mi){\ - dMiNi)(l - dLiM x N 2 ) 

As explained in the previous subsection, the generating function for products 

[k - Ai - A 2 , Ai, A 2 ] x[k- hi- V2, Mi, 0] x [k - fx 2 , ^ 2 , 0] (A.12) 
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is given by 



F 2 (d,L u L2,Mi,M2,N u N 2 ) = 

z R 

nnF^z' 1 ^ Li, L 2 , Mi, flf 1 , #2, M 2 , iVi, iV 2 ). 



(A.13) 



Here the variable z is introduced in order to keep the level fixed in the composition. By 
the determinantal formula, the generating function is essentially 



except that the coefficient of Mf 1 M^ 2 is the multiplicity of {jii — /j, 2 , jJ, 2 ). Thus the final 
generating function is 



This reproduces the generating function given in [4] and re-derived above. 

A. 3. Determinantal formula methods applied to the sp(4) andsu(4) cases 

In principle, the above procedure can be used to calculate the fusion rule generating 
functions for sit (4) and sp(4). Unfortunately, the intermediate expressions are too large to 
be manageable, even when manipulated with computer assistance. However, it is possible 
to calculate the specialisation of these generating functions with all but one variable, the 
level-grading variable, set equal to 1. For example, in the above calculation for su(3) we 
could have set L\ = L 2 = Ni = N 2 = 1 at the start of the calculation since they are not 
needed at any intermediate steps. Similarly we can set M 2 = M^ 1 at the last step which 
has the effect of setting M 2 = 1 in the final generating function. If we set all variables 
equal to 1, except the one that keeps track of the level, then the resulting generating 
function G(d) counts the number of independent couplings at each level. The su(4) and 
sp(4) specialised generating functions have been calculated in this way and the results are: 



F 3 (d,L u L 2 ,M u M 2 ,N u N 2 )= Q (F 2 -M 2 M^ 1 F 2 ) 



(A.14) 



F 4 = O F 3 (d, Li, L 2 , Mi, M 2 M^ 1 , Ni, N 2 ) 



(A.15) 



G su{4 \d) 



d 6 + 4d 5 + 13d 4 + 16d 3 + 13d 2 
(l-d) 12 (l-d 2 ) 



+ 4d+l 



(A.16) 



and 



G spW (d) 



d 4 + 2d 3 + 5d 2 + 2d+l 
(l-d) 9 (l + d) 



(A.17) 
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These expressions agree with the specialisation of the generating functions found in sections 
6 and 7; this thus provides a very strong independent verification of these results. In 
particular, it corroborates the closure of our set of fusion elementary couplings. 

Although we will not present the details of this derivation, we would like to draw 
attention to some technical issues. There are potentially two problems which could arise 
in using the deter minantal expansions. The first problem is that the determinant may 
contain terms which have level higher than the initial representation. For example in 
su(3) at level 1 the determinantal expansion of the representation (0, 1) is 



The representation (2, 0) is integrable only at level 2 and greater. However it can be shown, 
using the modification rules of [24] , that all such terms in the determinant vanish identically 
in the sp(2n) and su(n) fusion rings. Thus, when computing with the determinantal 
expansions at a given level, we need only consider terms corresponding to representations 
which exist at that level. 

The second complication which can arise is in a sense the converse of the first. There 
are representations which occur only at levels strictly greater than k, but which have 
determinantal expansions which contain products which are defined at level k. This does 
not occur for the su(n) determinants. However for sp(4) this problem can happen. The 
determinant formula for sp(4) is 



Take for instance the representation (0,2) which does not exist for level 1. However the 
determinant formula yields 



(0, 2) = (2, 0) ® (2, 0) - (3, 0) <g> (1, 0) - (1, 0) <g> (1, 0) + (2, 0) <g> (0, 0) (A.20) 



The only product which is defined at level 1 is (1,0) <E> (1,0) = (0,1). Thus the above 
determinant yields the following modification rule: (0, 2) = —(0, 1) for sp(4) at level 1 (see 
[24] for more details). Therefore, before converting the exponent of Mi into a Dynkin label, 
we must ensure that it is less than or equal to the exponent of d. This can be achieved by 
replacing M\ by Miy~ l and d by dy and then projecting onto non-negative powers of y 
and finally setting y = 1. 




(A.18) 




(A.19) 



39 



A. 4- Duality 



As described in [24] and references therein, there is a duality between fusion rules for 
su(n) at level k and fusion rules for su(k) at level n. This duality is somewhat involved 
when using standard Young tableaux. However, it can be clearly seen using contravariant 
tableaux. This duality can be used to provide a very nice nontrivial check of the sit (4) 
generating function. 

As discussed above, if all the grading variables in the su(4) fusion generating function 
are set equal to 1, except for the one associated to the level, we obtain (A. 16). However in 
order to use a duality argument to compare this expression with other generating functions, 
it needs some modifications. Duality maps Young tableaux to conjugate Young tableaux. 
For example su(3) at level 4 has 

(A.21) 



as a possible tableau and this maps to 



(A.22) 



in su(4) at level 3. In other words we have to include in the generating functions the terms 
corresponding to tableaux which have columns of length n in su(n). If J-'nid, Li, ...) stands 
for the original su(n) fusion generating function, then the procedure for incorporating 
tableaux augmented by columns of length n - while maintaining the first row smaller or 
equal to k — amounts to calculate 



9n(d) = 



d 2 



xyT n {dxy,x 1 L 1 ,x 1 L 2 ,...,y 1 M 1 ,y 1 M 2 , AT l5 N 2 , 



dxdy 

The effect of this operation is to multiply 



\x=y=l' 

(A.23) 



d k L?L? ...M^M^ ...N^N^ by (k - Ai - A 2 . . . + l)(k - ^ - H2- + 1) (A.24) 
which is the factor needed to add in all the Young tabl eaux with all allowed numbers of 



columns of length n. In other words, the su(3) tableau 
following equivalent forms: 



at level 5 should appear in 



(A.25) 
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that is, it should be counted four times. Doing this and setting all Dynkin-grading variables 
equal to 1 leads to the following generating functions: 

1 1 + d l + 3d + d 2 d 4 + 6d 3 + 10d 2 + 6d + 1 



90 = i 7r 91 = 7j 1^ 92 = — #3 = 

94 = 



1-d yi (1-d) 3 y2 (1-d) 6 yd ~~ (1-d) 10 
d 10 + 13d 9 + 78d 8 + 257d 7 + 513d 6 + 642d 5 + 513d 4 + 257d 3 + 78d 2 + 13d + 1 



(1-d) 12 (1-d 2 ) 3 

(A.26) 

The first two functions above correspond to the limiting algebras su(0) and su(l). For 
sm(0), there is only the trivial representation and it occurs at any level. Therefore, there 
is a single coupling at every level and there are no correction factors: go(d) = J2k d k ■ The 
function g\ can be constructed by duality. We start with the generating function for su(k) 
fusions at level 1. At level 1, we can ignore all relations between the elementary couplings; 
moreover, we can keep track only of those elementary couplings that occur at level 1: these 
are the various products involving the fundamental and the scalar representations. The 
truncated generating function then reads 

(A.27) 



(1 " d) nj(l - dLiNi)(l ~ dMiNi)} Y[ t j (l - dL^N^) 

where in the last series of term, the summation is defined modulo k with the understanding 
that Nk = 1. In this function, we replace d — > dxy, Li — > Li/x, Mi — > Mi/y, multiply the 
result by xy, differentiate with respect to x,y,d and set x = y = Li = Mi = Ni = 1, d = 
(to keep only the linear term in d). This gives (k + l) 2 . Hence we have 

9i(d) = J2( k + 1 ) 2 d k = 7 T ^s (A.28) 

These functions g n (d) display very nice properties: 

1- the factor (1 — d) occurs to the power (n + 2)(n + l)/2 in the denominator; 

2- the numerator polynomial p n (d) satisfies p n (l/d)d deg ( pn ^ = p n (d); 

3- p n (d) has positive coefficients; 

4- the difference between the degree of the numerator and denominator is In. 

The mere fact that g±{d) shares the generic properties of the previous g n functions is 
supporting evidence for the correctness of the sit (4) generating function. 
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The Taylor expansions of the g n (d) functions read: 

g (d)= 1 + d + d 2 + d 3 + d 4 + d 5 +■■■ 
gi(d)= 1 + Ad + 9d 2 + 16d 3 + 25d 4 + 36d 5 +■■■ 
g 2 (d)= 1 + 9d + 40d 2 + 125d 3 + 315d 4 + 686d 5 +••• 
g 3 (d)= 1 + 16d + 125d 2 + 656d 3 + 2646rf 4 + 8832rf 5 + ••• 
g 4 (d)= 1 + 25d + 315d 2 + 2646d 3 + 16720rf 4 + 85212c/ 5 +••• 

(A.29) 

from which duality (i.e. horizontal versus vertical) is completely manifest. (We stress 
that the 'built-in duality' for obtaining gi concerns only the second row and the second 
column.) In particular the first 4 terms 1, 25, 315 and 2646 of the su(4) function match 
the coefficients of the 5-th column. This again provides independent evidence for the 
correctness of the su(4) fusion generating function out of which the function g^ has been 
constructed. In particular, this is a decisive test of the necessity of the extra elementary 
coupling F and an evidence for the absence of further additional elementary couplings. 

From the above functions g n (d) we can construct the sum 

f(r, d) = g (d) + 9l (d)r + g 2 (d)r 2 + g 3 (d)r 3 + ... (A.30) 

where r is the grading variable associated to the rank +1 (i.e., its exponent is the value of 
n for su(n)). It satisfies f(r,d) = f(d,r) by duality. We speculate that other symmetry 
properties might be used to provide an explicit formula for /(r, d). 

We can illustrate this dual symmetry in a particular example. Consider the function 
g n (d) that counts the number of couplings of the representation [k — 1, 1, 0, 0, . . . , 0] with 
anything in su(n) at level k. Since the Young tableau of (1, 0, 0, . . . , 0) is invariant under 
a duality transformation exchanging k and n, by summing up the resulting functions 
multiplied by r n , one should produce an expression f(r,d) symmetric in the interchange 
of r and d. The function g n (d) is calculated as follows in terms of the original su(n) fusion 
generating function T n : 

d 2 

g n {d) = q m ^q x xF n {dx,x 1 L 1 ,x 1 L 2 , Mi, 1, 1, l...)U=y=i,Mi=o,i,i=...=i (A.31) 

As explained above, the differentiation with respect to x is required in order to take 
into account all contributing Young diagrams associated to the first representation (A). 
The second representation being fixed to be (1, 0, 0, . . . , 0), does not require an adjusting 
multiplication factor. Setting the variable Mi = 0, after having differentiated with respect 
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to it, simply serves to select the term linear in Mi. Since the representation (1, 0, 0, . . . , 0) 
does not exist for s-u(O), go(d) = 0. The function gi(d) is found by duality as explained 
previously. The first few su(n) functions g n (d) are found to be: 

2d - d 2 ~ / r\ 3d - d 2 

_ , A U - d 2 _ . „ 5d - d 2 v ' 

9M = jj-^ 94 (d) = j^—^ 

Fortunately, the general pattern is clear: the expression of g n is easily guessed to be: 

(n + l)d- d 2 

9n(d) = (1 _ d)n+1 n > 1 (A.33) 
From this exact form of g n (d), we can write down readily the exact expression for the sum 

/>, d) = ± Ud)r- = d £ 2 ~ d _~$ (A-34) 

n=l { ' 

The result is manifestly invariant under the duality transformation that interchanges r and 
d. 



Appendix B. Status of previous conjectures 

In this appendix, we would like to clarify the relation between the present work and 
our previous ones and state precisely in what sense our previous conjectures are either 
embodied in the present reformulation of the problem or have been proved. 

A general approach to the construction of generating functions for fusion rules was 
proposed in [4]. It was based on the following two conjectures: 

1) Every coupling is characterised by a threshold level k . The multiplicity of a triple 
product at level k is given by the number of couplings with threshold levels < k. 

2) There is a choice of forbidden couplings such that the threshold level of a coupling is 
the sum of the threshold levels of its components. 

As already mentioned, it can be shown [6] that conjecture 1 is a consequence a sharp- 
ened formulation of the depth rule [26]. This leaves us with a single conjecture which we 
rename: 
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Conjecture I: There is a choice of forbidden couplings such that the threshold level of a 
coupling is obtained from the sum of the threshold levels of the elementary couplings that 
appear in its decomposition. 

In the formulation of conjecture I, the element of 'choice' refers to the fact that both 
sides of a tensor- pro duct relation do not always have the same threshold level and which 
one is taken as the forbidden coupling makes a difference in the generating function for 
fusion rules. With the notion of a set of elementary fusion couplings, which includes the 
scalar one (this is a new feature of the present work), all relations acquire equal threshold 
levels and this choice becomes immaterial. This suggests the following modification of 
conjecture I: 

Conjecture V: The threshold level of a fusion coupling is read off from its decomposition 
into the elementary fusion couplings. 

A interesting aspect of this reformulation of the conjecture is that it embodies an 
observation that was presented as a conjecture in [7], namely that the level is always 
minimised. More precisely, in the choice of forbidden couplings, we should always forbid 
the one with higher threshold level. This 'minimal level' prescription is automatically 
taken into consideration here since the relations have identical levels. If one of the products 
appears in the relation with a factor E , it means that the product without this E factor 
occurs at a lower level and it is not forbidden. For instance, the relation E1E3E5 = 
EoEiEs indicates that the coupling E^Eg appears at level 2. In the tensor-product relation 
EiE 3 E 5 = E^Eg, we have thus effectively forbid the higher-level term of the relation. 

Once the notion of fusion elementary couplings in terms of which every coupling can 
be decomposed (conjecture I') is introduced, this naturally calls for a reinterpretation in 
terms of a fusion basis. It is indeed plain that our conjecture (and the mere existence of 
threshold level) boils down the fundamental conjecture presented in the text, that is, the 
existence of a fusion basis. 
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